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1 Introduction 



Of late there has been a great deal of interest in discrete integrable systems, owing to 
their numerous interesting properties. One such property which has emerged lately, is 
the existence of the so called spectral curve associated with an integrable system[l]. 
The latter in turn is closely connected with the motion of separability. Indeed, it 
is well known that the technique of separation of variables provides an invaluable 
tool for the analytic construction of action-angle variables of an integrable system [2]. 
Furthermore, it may be shown that the locus of points, given by the separation 
variables, define in fact the spectral curve. 

On the other hand, the concept of an r-matrix algebra marks a major development 
in theory of nonlinear integrable systems especially in regard to their Hamiltonian 
nature [3]. It provides an elegant formalism for proving the Poisson involutiveness 
of the integral of motions for such systems. However, it was only in the eighties 
that the connection betweenthe r-matrix formalism and separation variables of such 
systems clearly emerged as a result of the efforts of Komarov [4], Sklyanin [5] and 
other researchers of the St. Petersburg School. 

In this communication, we examine the nature of the relationship between the 
separation variables of a discrete lattice integrable system and the corresponding 
local lattice variables. In general, the nature of the transformation relations between 
these two sets of variables is rather complicated and an explicit determination of 
the generating function proves to be rather daunting. The existence of the latter is 
essential to gaurantee the canonical nature of the transformation. However, recently 
Kuznetsov [6] has shown that, not only does the inverse transformation between the 
separation variables and the local variables exists, but also that these trasformations 
can even be broken down into a sequence of transformations each with a generating 
function. 

The organisation of the paper is as follows. In section 2, we describe the basic 
features of the discrete model under consideration. The r-matrix algebra of the mon- 
odromy matrix is explicitly solved, as also the leading order structures of its elements. 
In section 3, we introduce the spectral curve and separation representation for the 
elements of the monodromy matrix. In section 4, which constitutes the main body 
of the paper, we derive explicitly the inverse separating map between the separation 
variables and the local lattice variables defining our system. The canonical nature of 
this mapping is displayed by working out the corresponding generating function. 

2 Formulation 

We begin by considering a discrete integrable system which is described by a local 
Lax operator l n (u) of the form 




(2.1) 
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Here u is the complex spectral parameter and q n and p„are the nonlinear variables 
defined at the nth lattice points. They satisfy the canonical Poisson brackets 



{Qk,qj} = {Pk,Pj} = for all k,j = 1,2, ,n and {qk,Pj} = $ kj 

It is a matter of simple computation to show that 

{l n (u) ® l m (v)} = [r(u - v),l n (u) <g> l m {v)\8 mn 
where r(u — v) represents the classical r-matrix given by 



(2.2) 
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with P being the permutation matrix satisfying P 2 = I. We define the monodromy 
matrix T?q{u) as the ordered product of l n {u)'s over the entire lattice. That is, 



T N (u) = J| l n (u) = In^In-^u) l2(u)k(u) 

n=N 

Upon using ( 2.3 ) it is easy to deduce that 

{T N {u) f T N {v)} = [r(u - v),T N {u) <g> T N (v)} 



(2.4) 



(2.5) 



Consequently we infer that structurally the Poisson algebra of l n (u) and Tjv(u) are 
identical. In the literature on nonlinear integrable systems equation (2.5) is usually 
referred to as the quadratic r- matrix algebra. In case of 2 x 2 local Lax operators 
the monodromy matrix is itself a 2 x 2 matrix which may formally be written as, 



t n{u) =H l n (u) 



n=N 



A N (u) B N (u) 
C N (u) D N (u) 



(2.6). 



The Sklyanin algebra ( 2.5) then takes the following explicit form, namely 

{A N {u), A N (v)} = {B N (u), B N {v)} = {C N (u),C N (v)} = {D N {u), D N {v)} = 0, 



{C N (u),A N (v)} 
{A N {u),D N {v)} 
{B N (u),A N (v)} 
{C N (u),D N (v)} 
{B N (u),D N (v)} 
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{A N (u)C N (v 
{C N {u)B N {v 
{B N {u)A N {v 
{C N {u)D N {v 
{D N (u)B N (v 



A N (v)C N (u)}, 
C N (v)B N (u)}, 
B N {v)A N {u)}, 
C N (v)D N (u)}, 
D N {v)B N {u)}. 



(2.7) 



From the definition of monodromy operator in (2.6) and the form of the local Lax 
operator as given by (2.1) by explicit calculation obtain the leading order terms of 
the matrix elements An(u), B^{u) etc . These are as follows 

a n {u) = u N + e:^<i,i>,)u x 1 + 

(2.8a) 

B N {u) = ( ^ =lPi )u N ; + 

(2.86) 

Cn(u) = (El^pju"- 1 + 

(2.8c) 

D N (u) = u N -(T l f =1 q i p i )u N ' 1 + 

(2.8d) 

Note that detl n {u) = u 2 so that detT N {u) = u 2N . Now if at each lattice site an 
inhomogenity parameter c n is introduced, then the local Lax operator assumes the 
following form: 

and consequently the monodromy matrix 

m ( f \\ I A n(u) B n (u) \ 
T»(«,{c»})=[ Cn{u) Dn(u) ) 

now involves the inhomogenity parameters {c n }^ =1 . Its elements then are as given 



(2.10a) 
(2.106) 
(2.108c) 
(2.10d). 

Thus in general A^,Dn are polynomials of degree N, while .B^and Cn are polyno- 
mials of degree N-l in the spectral parameter. Furthermore, in this case 

det l n {u - c n ) = (u - c n ) 2 

so that 

N 

det T N (u, {c n }) = J] (u - c n ) 2 (2.11) 

71=1 

The involutive, independent integrals of the system are now obtained as the coeffi- 
cients of trace of the monodromy matrix. That is, 

trT N (u, {c„}) = A N (u, {c n }) + D N (u, {c„}) 



below 






A N (u, {c n }) 


= u iV +(S£ 1 g i p i -E£ lQ 


)^- x + 


B N (u, {c„}) 


= (-i:l 1 p l )u N - 1 + 




C N (u, {c n }) 


= (^ =iq ^)u N - 1 + 




D N (u, {c n }) 


= u N -(T,* =1 q i p i +Y 1 *L 1 Ci 


)^+ 
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= 2u N + Hm"- 1 + H 2 u N ~ 2 + + H N _ lU + H N (2.12) 

For a three particle system i.e N=3 

E x = -2( Cl +C2+c 3 ) 

H 2 = 4:(c 3 c 1 + 1 C2+C2C 3 )+2(k 3 k2+k 2 k 1 +k 1 k 3 )-k 3 p 3 p2-k 3 p 3 p 1 -k2P2P3-hP2Pi-kiPiP3-kiPiP2 

Note that throughout this paper we shall assume periodic boundary conditions i.e. 

q n+N = q n and p n+N = p n 

3 Spectral curve and Separation representation 

In case of Liouville integrable systems, the technique of separation of variables play a 
vital role in construction of analytic expressions for the action angle variables. How- 
ever it was only in the early eighties that the connection between separation variables 
and the r- matrix formalism of the quantum inverse scattering method became evident. 
The fact that for 2 x 2 Lax matrices, the separation variables are given by the zeros 
of the off diagonal of T^{u) was first observed by Komarov [4]. The technique was 
investigated in great detail by Sklyanin [5] and was applied to several models, such 
as the Goryachev-Chaplygin top, the Neuman system [8] etc. Indeed, its subsequent 
application to the 0(4) Kowalevski top, by Kuznetsov, represents the culmination in 
a sense, of this remarkable technique[9] . 

We shall briefly outline the connection between separation variables and the el- 
ements of the so called spectral curve, which is explained below. It is well known 
that, an integrable Hamiltonian system possessing N independent integrals of motion 
in involution, is separable, if there exists a set of Darboux co-ordinates say (itj, Vi)^ =1 
which satisfy a set of equations of the form 

f{m,Vi,h,I 2 , ,In)=0 (i=l,2, ,iV) (3.1). 

These equations arise from the equation of an associated spectral curve T defined by 
the following equation 

det(zl - T N (u)) = z 2 - P(u)z + Q(u) = 0, (3.2) 

on the ( u, z ) plane. Here P{u) and Q(u) are given by tr(Tjv(it)) and det(TN{u)) 
respectively. In our case we note that Bn(u) i.e the upper off diagonal element of 
Tjv(m) is a polynomial of degree (N-1), so that it can have at most ( N-1) distinct 

zeros, which we shall denote by u, L (i = 1,2, , N — 1). One may show that there 

exists a system of differential equations for these zeros. It is easy to see that at the 
zeros of B N (u), viz. Ui (i = 1, 2, N— 1) the monodromy matrix becomes triangular, 
i.e, 

rp ( \ ( A N (Ui) \ , 

Tn{Ui) = { C N {u t ) D N (ui) ) ^ 
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and therefore its eigen-values are the diagonal elements A N {ui) and D N {ui) repec- 

tively. If Zi = A N {ui), i — 1,2, , N — 1, then it is obvious that the set of points 

(ui, Zi)f =l of the (u-z) plane must lie on the spectral curve Tj 

zf - P(u t )z, t + Q( Ui ) = 0^ = 1,2, , N - 1. (3.4) 

For the remaining part of this communication we shall set Zj = e _1,i and shall show 
that the pair (ttj, Uj) constitute a set of Darboux co-ordinates satisfying the canonical 
brackets 

{u i ,u j } = {v i ,v j } = 0, {ui,Vj} = S i:j i,j = 1,2, ,7V-1. (3.5) 

Let us further define two quantities as follows: 



u n = ' v ^ = -l. i=1 Pi, 



satisfying {un, v n} = 1 Thus altogether we have a total of N pairs of Darboux co- 
ordinates {ui,Vi)f =1 . Having defined the above separation variables we next proceed 
to the construction of a representation of the quadratic r-matrix algebra (2.5) in terms 
of these separation variables. This is easily accomplished by using Lagrange interpo- 
lation. For instance, we may express B N {u) in terms of its zeros u±, u 2 , un-i 

as 

B N (u) =v N Y[(u-Ui). (3.7) 

n=l 

Furthermore, assuming these zeros to be mutually distinct we may express any poly- 
nomial F(u) of degree less than (N-l) in the following mannner 

FM Sgi^ -BTT^ ^ (3-8) 



B N (u) ' B' N (ui)(u - UiY 
while if F(u) has degree > (N-l), then 



F{U) = S£r 1} n , ; x + G(u) (3.9) 



B N (u) 1 1 B' N (ui)(u - Ui) 

where G(u) is a polynomial of appropriate degree. In the present case since An(u) 
and Dn(u) are polynomials of degree N, we have the following representations for 
them, 

A N {u) - B N {u)[ + S l=1 B , N{ui){u _ Ui) \ (3-10) 

D N (u) = BN{u)[ --(-NV N -^ 1) u i + 2i;l 1 c t ) + detT N ( Ui y< 

v N B N {ui){u- Ui) 

since 

detT N { Ul ) = A N ( Ui )D N ( Ui ) = e- Vi D N ( Ui ) (3.12) 
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The representation for Cn(u) follows from the relation 

detT N (u) = A N (u)D N (u) — B N (u)C N (u) 

so that 

r / x A N (u)D N (u) - detT N (u) 

Cn{u) ~ (3 - 13) 

with An(u), Dn(u), and Bn(u) as given by the relations (3.10) , ( 3.11), and (3.7) 
while 

1 N 

detT N = Yl detl n (u - c n ) = J| [u - c n ) 2 . 

n=N n=l 

We shall now proceed to show that the above representations for the elements of 
T N {u) does indeed reproduce the quadratic r-matrix algebra. For example, 

{B N (u), A N (v)} = {B N (u),B N (v)[ V + {UNVN + ^ 1)u ^ + Sf--— ,]} 



v N 1 B' N (ui)(v -m) 

B N (v){B K („), » + («^ + ggr"«.,) + -} 

v N B' N ( Ui )(v - u^ 1 



since {B N (u), B N (v)} = 0. 



N-l i=l -y N-l 

= B N (v){v N ,u N } Yl (u-Ui)+B N (v)T,^ =1 -lv N Y[(u-u k ){-Ui,e v '} JJ (u-u k ) 

n=i k=i t> N [Ui)[y - Ui) k=i+1 

B N {v)B N {u) B N {v)B N {u) N _ x e-* . 1 1 . 

H ~ — ~ ui („. \ ^7, — ~ — 7. — 77> 



vn u — v 1 B' N (ui) v — Ui u — Ui 

B N (u)B N (v) B N (u)B N (v) A N (v) v + u N v N + Y^i^Uj A N (u) u + u N v N + E^ -1 ^) . 



vn u — v Bn(v) vn Bn(u) vn 

{B n (u),A n {v)} = -±-(B N (u)A N (v) - B N (v)A N (u)). 

It is also straight forward to show that detTN(u) is a Casimir, i.e, its Poisson brackets 
with the elements of Tn(u) vanish. 



4 Factorized separation 

From the above discussion it is obvious that we have at our disposal two sets of local 
variables- the first set consisting of (qi,Pi)iLi, in terms of which the local Lax matrices 
are represented and a second set consisting of the separation variables (ui, Vi)f =l . The 
two sets of variables being connected by the separation map Si given by 

B N (ui) = 0, e Vi = A N (ui) i = 1,2 ,(N-1) (4.1a) 
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while 

wat = —p , v N = -E i=1 pi. (4.16) 

The concept of factorized separation has recently been introduced by Kuznetsov in 
[6] to obtain an inverse of the separation map. Indeed, the mapping Sn between 
the original local variables ijpi,qi) and the separation variables is in effect a rather 
complicated canonical tranformation and to obtain its generating function is a tedious 
task. The same also holds for the inverse mapping which relates the separation 
variables to the local variables appearing in the Lax operators. In [6] Kuznetsov 
made the remarkable observation that not only does the inverse mapping S^ 1 exist 
but that it could also be factorized as a chain of more simpler, at least in principle, 
canonical transforms, thus leading to the solution of the inverse problem. 
We note that the N-site monodromy matrix may also be expressed as 

T N (u) = l N (u - c N )T N ^(u) = ( u - c " + 9nPn -P» ) T^u), (4.2) 

^ q%p N u-c N - q N p N J 

while 

Tat_i(m) = l^(u - c N )T N (u) = /j^f U n , T N (u). (4.3) 

aetlN{u — cn) 

If we formally denote 

then it is obvious that Bn-i{u) is a plynomial of degree (N-2) while An-i(u) and 
Dn-i(u) are polynomials of degree (N-l) in u. Moreover, by analogy with their earlier 
counterparts, let Uj (j — 1, 2, , (N — 2)) denote the zeros of -Bjv-i, so that 

B N - 1 (u)(u j )=0, j = 1,2, ,(N-2) (4.5) 

and let 

e~ Vi = A^fa); j = 1,2, ,(N — 2) (4.6) 

while 

Sj =1 (qiPi — Ci) _ v,AT-l (a ^7^ 

un-i = ^iv^T , v n-i = -S i= i Pi (4.7) 

~^i=l Pi 

As before A N _i(u) and D N _ 1 {u) may in view of (4.5) and (4.6) be expressed by 
Lagrange interpolation formulae in the undermentioned forms: 

/ x / x < u + {u N v N + E^r^tL) , r „ e~ Vj 

A K -M = Bs-Ml + ^f ^— -y l (4.8.) 

L>n-i\u) = i^jv-iW ~ 1- E =1 — — — — 

(4.86) 
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while 

JV-l 

B N -i(u) = vn-i H(u- Uj). (4.9) 
j'=i 

Let Xn denote the mapping from the separation variables itj, = 1,2, , (TV — 1)) 

to the new separation variables Uj,Vj (j = 1,2, , (TV — 2) which parameterize 

Tn-i(u) and into the pair of local variables and ^tv- From (4.3), since det In(u — 
cat) = (m — c^) 2 , we see that Tjv-i(m) = ^^z^t^ ^V (~») has a second order pole at 
u = cat. As Tjv-i(m) does not have any singularities, consequently upon equating the 
residue at u = c^- to zero, we have 

Jrn^ ^-[Adj l N (u - c^TVH] = 0, (4.10) 
which leads to the following relations: 

( -QnPn Pn \ = L( A n(cn) B n {c n ) \ ( £D N (u) -£B n {u) \ 
V -Q 2 nPn QnPn J A \ C N {c N ) D N (c N ) J \ -£C N {u) £A N (u) 

where 



(4.11) 



Hence 



and 



Pn = ^[-A N (c N )^j-B N (u) + B N (c N )^-A N (u)] u=CN (4.13) 

„ r C N (u)£B N (u) + D N (u)£A N (u) 

qN [ -A N {u)±B N {u) + B N {u)±A N {uY u=CN l4 - 14J 

Equations (4.13 - 4.14) determine the Nth local variables in terms of the separation 
variables (itj, Vi)\f =v To determine the remaining local variables we adopt a recursive 
procedure. To this end we shall first consider the mapping 

X N :(u,v)\ — > (u,v\q N ,p N ), (4.15) 

where u — (ui, 112, ujv)eic. Next if we consider the mapping X^-i given by, 

Xat_i : (u,v\q N ,p N ) 1 — > (u,v\qN,PN;qN-i,PN-i) (4.16) 
then it is obvious that the composition of these two mappings, 

Xn^oXn) : (u,v) 1 — > {u,v\q N ,p N ;q N ^p N ^) (4.17) 
Therefore the composition of N such mappings viz 

X l0 X 2 o o X N _! o X N : (u, v) 1 — > (q,p) (4.18) 



8 



leads to a transformation from the separation variables {ui,Vi)\f =l to the local vari- 
ables (qi,Pi)iLi , which is precisely the definition of S^i.e, 

S N 1 = X 1 oX 2 o oX N -!oX N (4.19) 

Equation ( 4.19) represents the separation factorization chain of the inverse mapping 
from the separation variables to the local variables occuring in the Lax operators. 

In the present case the mapping X N will be explicitly derived in the following 
manner [6]. From equation (4.2) we obtain 

B N (u) = (u- c N + q N p N )B N ^i(u) - p N D N _i(u) (4.20) 

If this equation is evaluated at u — — 1,2, ,N — 2) i.e, at the zeros of B N _i(u) 

then one obtains the following relation 



(4.21) 
(4.22) 



B N (uj) = -PnDn^uj) j = 1, 2, , AT-2 

B N (uj) = -p N {detT N ^(u 3 )e^} j = 1,2, ,N-2 

N-l 

= -pn n - c kf^ 3 

k=l 

Again from (4.3) we obtain 

B N -i(u) = / r[(u - c N - q N p N )B N (u) +p N D N (u)] (4.23) 

(let In(u — cjv) 

Its evaluation at u = Ui gives 

{u t -c N fB N ^ l {u i )=p N D N {u i ) 2 = 1,2, ,N-1 (4.24) 

Substituting expressions for B N {u) and Bn-i(u) from (3.7), and (4.9) respetively, 
and D N -i(u) from (4.8b) into (4.20) and equating the leading order terms, we obtain 
from the coefficient of m^" 1 , 

v N = v N ^! -p N . (4.25) 

However, it is also follows from the definitions of (un,vn) and (un-i,vn-i) as given 
in (4.1b) and (4.7) repectively that 

u N v N = T^ =l (qiPi - a) and u N -iv N -! = E^Li^qy?* - q), (4.26) 

so that their difference i.e, 

U N V N - Ujv-l^n-l = qNPN - C N . (4.27) 

Now from (4.25) and (4.27) we may solve for vn and vn to obtain, 

PnU N -i + q N p N - c N 
v N = (4.28a) 
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PnU N + q N p N - c N 
vn-i = (4.286) 

UN — V-N-i 

Equating the coefficients of u N ~ 2 in (4.20) we obtain 

-VatS^ 1 = VN^i-Zfj^Uj + (q N p N - c N )} +p N (u N v N + 2E^ 1 c i ) 
leading to 

-V N Z i=1 Ui + V N -i{H j=1 Uj + C N ) 

{u N v N + 2E i=1 a + q N VN-i) 

Substituting the expressions for vn and vn-i from (4.28a,b) into (4.29) yeilds the 
following quadratic expressions for Pn(qn, u, u, {cj}) namely, 

P 2 n (Qn + u N )(q N + u N -!) +p N [(qN + UN-\){^iJiUi - c N ) - {q N + u N )(T,fs 2 Uj + c N ) 

+2Eg- 1 1 c i (u N - Uiv_i)] + c N [c N + ?%-?Uj - EjllV] = 0. (4.30) 

Incidentally, if we set q = for all i — 1, 2, AT, then (4.30) leads to the following 
linear expressions for 

pn = [ ^ - ( 4 - 31 ) 

gAr + u N _! q N + u N 
On the other hand from (4.22) we get 

-v 3 _ Pn UkJi (Uj-c k ) 2 _ p N Uk=i (uj-c k ) 2 _ , 

Upon substituting from (4.28a ) we get 



PnKQN + UN-i) - C N ^ [Uj ~ U k ) 

Finally from (4.24) we have upon using the expressions for B N _i(ui) and D N {uj) 
obtained from (4.9) and (3.11) repectively, the following relation, 



PN I\m=\{ u i~ c mf 
V N -l(Ui - Cat) 2 nfj 2 (Ui - Uj) 



e-* = , Z : - 2 , i = l,2,...,JV-l (4.34) 



whence using (4.28 b ) to eliminate v^-i we get 

e -vi = Pn(u n - Mjy-ij 1 YlZzljuj -c m f i = 1 2 N — 1 

PN{qN +«JV-l) - Cat (Wj-Cat) 2 nft 2 («i ~ «j) ' 

Thus the mapping Xat is formally defined by the following relations 

Pn(Qn + u N )(q N + Ujv-i) +p N [(qN + UN-i)(^'i f =i 1 Ui - Cat) - (q N + u N )(Hf~ 2 Uj + c N ) 
+2Z? = ~ 1 1 c l (u N - ujv-i)] + c N [c N + Tff~ 2 Uj - T^Ui] = (4.36) 



(4.35) 
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e - n= V»W N -u N1 ) l^ gjgfe^g i = 1 , 2 ,.... iJV _ 1( 4.37) 

PM<?JV + UN) - C N [Ui - C N y Uj=l(Ui - Uj) 

and 

e - h = _ Mu„-u N1 ) "tffp^ctf N _ 2 

Pn(,QN + - Cat ^ [Uj - u k ) 

together with the relation 

Pn{u N -i + Qn) - c N _ Pat(mat + gjv) - , . 

vtv = = , ujv-i = = (4.39) 

U N - «n-l WAT - U n _i 

For the special case when q = 0, for alii = 1,2, , JV these expressions reduce to 

the following 

^f=i 2 Uj ^fjiUi Pn(u N -i + gjv) _ Pn(u n + Qn) 

Pn = [ — ; , v N = - , v N -i = - 

q N + u N ^i q N + u N u N - m A t_ 1 u N - u N _ x 

(4.40a, 6, c) 

e-* = " ^ < = 1,2, ,iV-l. (4.40d) 

e~*> = - UN -y- 1 UjJ^—, J = 1,2, ,N-2. (4.40e) 

OAT + Mat_i ^ (Mj - U fc ) 

It is easy to verify that the expressions given in equations (4.40 a,b,c,d,e ) are derivable 
from the following generating function F(qx] u, u), by means of the following mapping 
X N : 

dF 

Vi = Tr- i = 1,2, , N. 

OUi 

dF 

*' = -^ •< , n — i. 

dF 

Pn = -j—. (4.41) 
dq N 

with 

F(q N ; u,u) = [-H^Uj]log(q N + u N -i) + [E^VlW^iv + u N ) 



+ [S, =1 - S i=1 Ui}log(u N - 



+2(iV - l)^Mu 3 ) - 2(iV - 2)Sf = I 1 6(«,) 

+E^- 1 Ef = - 1 2 e(^ - uj-) + m(N - 2)Ef = - 2 ^ (4.42) 

where O(m) = / log{u)du + c(constant) . The generating function for the subsequent 
mapping .Xjv_i is obtained by replacing N by N — 1 and so on. The above formulae 
thus describe explicitly the inverse separation map S^ 1 and by virtue of the existence 
of a generating function for each stage we conclude that the inverse separation map 
is canonical. 
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5 Discussion 



In this communication we have studied the nature of the inverse separation transfor- 
mation between the local nonlinear variables of a discrete integrable system, and the 
corresponding separation variables of the system. The novel feature of the method- 
ology used is the fact that this separation transformation can be factorized into a 
sequence of transformations, each with its corresponding generating function. The 
latter clearly is indicative of the canonical nature of the separation transformation. 



6 Acknowledgement 

One of the authors (S.M), is grateful to UGC for a Project Associate fellowship which 
made this work possible. 



12 



7 References 



1. V.B. Kuznetsov, M. Salerno, E.K.Sklyanin J.Phys. A. Math Gen.33A (2000), 171 

2. V. B. Kuznetsov, M.Petrero, O.Ragnisco, \nlin. SI/040302% (2004) 

3. L.Faddeev and L. Takhtajan, Hamiltonian methods in the theory of solitons, 
Springer Verlag, Berlin, (1988) 

4. I.V.Komarov, Theo. Mat. Fiz., 50 (3), 402, (1982). J.Phys.A15 (6) (1765) 1982 

5. E. K. Sklyanin Prog. Theo. Phys. Suppl. 118 , 35, (1995) 

6. V. B. Kuznetsov -Inverse problem for sl(2) lattices nlin.SI/ 0207025 

7. K.Takasaki, Spectral curve, Darboux coordinates and Hamiltonian structure of pe- 
riodic dressing chain, nlin.SI/ '02060^9 

8. E. K. Sklyanin, Differential Geometry, Lie groups and Mechanics, VI Zap. 
Nauchn. Sem. Leningrad, (LOMI) 133, 236, (1984). 

9. V. B. Kuznetsov, Kowalevski top revisited, SI/0110012. 



13 



